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ABSTRACT 

We present a numerical study of Gamma-Ray Burst - Circumburst Medium interaction and plasma precondi- 
tioning via Compton scattering. The simulation tool employed here - the PhotonPlasma code - is a unique 
hybrid model; it combines a highly parallelized (Vlasov) particle-in-cell approach with continuous weighting 
of particles and a sub-Debye Monte-Carlo binary particle interaction framework. Our first results from 3D sim- 
ulations with this new simulation tool suggest that magnetic fields and plasma density filaments are created in 
the Wakefield of prompt gamma-ray bursts, and that the photon flux density gradient has a significant impact on 
particle acceleration in the burst head and wakefield. We discuss some possible implications of the circumburst 
medium preconditioning for the trailing afterglow, and also discuss which additional processes will be needed 
to improve future studies within this unique and powerful simulation framework. 

Subject headings: Gamma rays: bursts - plasmas - magnetic fields - methods: numerical - instabilities - 
acceleration of particles 
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1. INTRODUCTION 



The microphysics of Gamma-Ray Burst (GRB) afterglow 
shocks propagating through a circumburst medium (CBM) 
has been theoretically predicted by (Medv edev & Loeblll999|) 
and later extensively studied using p article-in-cell (PIC) 
models in the past few years (see e.g. iKeshet et al.l 12008 



iHededal et al."l26o4j; iFrederiksen etani2004l: ISilva et al.ll2003 , 
Kjruzinov 
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In contrast, the initial interaction between 
the GRB photons and the CBM plasma is less well studied 
numerically, although progress has been made for coherent 
wake field processes in lower-dimensional an d pair plasmas 
(e.g. Hoshinol 120081: iLiang & Noguchill2007h . Whereas the 
burst-C BM interaction has r eceived strong attention in the- 
oretical (Beloborodov 2002; Thompson & Madau""2000) and 
observational ( Connaughton 2002) works, the detailed micro- 
physical plasma dynamics in the first pulse-plasma encounter 
lacks sufficient treatment with respect to the high degree of 
forcing anisotropy in the problem of photon scattering off the 
relatively cold tenuous CBM. 

Wakefield interactions between high energy photon pulses 
and a plasma have been extensively studied in the context 
of coherent wakefield ac celeration in laser-plasma systems, 
both experimentally (e.g. lPhuoc et al . 2005; Kostvukov et al. 
| 2003 | ) , and num erically in the same context (e.g. Phuoc et al. 
120051: IPukhoV 2000). In the astrophysical context of parti- 
cle acceleration in GRBs, wakefield acceleration has been 
studied in the coherent regime where the forcing of the 
plasma is done by the e lectromagnetic field (iHoshinol 120081: 
iLiang & Noguchill2007l) . Incoherent - or stochastic - wake- 
field processes have to our knowledge so far not been studied 
numerically, a nd has been considered the oretically only in a 
few cases (e.g. lBarbiellini et al.ll2006L 12004 ). 

Our motivation for developing a new and improved 
Particle-in-cell- Monte-Carlo (PIC-MC) code framework for 
investigating the detailed Compton interaction of an ultra- 
intense GRB and a quiescent cold CBM is two-fold: 

First, due to the high-energy spectrum of the prompt GRB 
(for BATSE of order ^ nieC^), the photon wavelength is 
much shorter than any characteristic length scale in the CBM 
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plasma, <C Se and <C A^, where Se and A^ are the elec- 
tron skin-depth and Debye-length, respectively. We there- 
fore expect a binary approximation to work well for prompt 
burst photons interacting with the plasma constituent parti- 
cles through Compton scattering. In fact, photons in a large 
range down towards the electron plasma frequency uj^ ~ ujpe = 
^r)^th,e (where Vth^e is the electron thermal speed) could be 
treated as point particles that interact in a binary way with the 
plasma. 

Secondly, the extreme intensity of the burst makes the GRB 
front opaque to the plasma electrons. Assume a CBM with 
ne ^ I cm~^ , and f urther assume a typical BATSE GRB with 
an estimated (from lBarbiellini et"aDl2006l eqn.2) photon flux 
energy density of 10^^ eV cm~^ s~^ at distance Rq ~ 10^^ cm 
from the progenitor. Then, Compton scattering off the CBM 
plasma will occur many times every second per electron dur- 
ing the prompt phase; most photons survive traversing the 
plasma without being significantly affected. The plasma is 
forced predominantly by binary particle- particle encounters, 
everywhere locally. 

These considerations motivated us to develop and employ 
a new and unique tool in order to gain access to the sub- 
Debye j;egime;__Adopting the "random phase approxima- 
tion" ( Pines & Bohmlll95 2) we may split, in a natural way, 
plasma processes into two different numerical schemes de- 
pending on the characteristic scale of the dynamics, here de- 
noted kc: 

Az) <C |kc|^: for forcing and disturbances scale larger than 
the Debye length the Vlasov approach is applied and 
the PIC code framework is utilized. 

A/) |kc| S in the opposite case dynamics is particle- 
particle interaction dominated, and we solve the scatter- 
ing problem statistically through MC radiative transfer 
for the plasma constituents. 

In this way we obtain a substantially improved description of 
the plasma for both the Vlasov and detailed balance cases, 
and capture in a natural way the details of wave-particle 
interaction and scattering. The PhotonPlasma code is 
parallelized with MPI and is highly scalable. Any plasma 
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constituent species, forcing and scattering process is possi- 
ble to incorporate with relative ease; for further details see 
(lHaugboellell2QQ5l:lFrederiksenll2QQ8h . 

In the remainder of this Letter we present, in Section 2, 
some central details on our numerical approach, and account 
for the simulation setup as well as physical and numerical 
scaling. Results and discussion are given in Section 3, where 
we also relate the wakefield effects to other results concern- 
ing coherent - rather than stochastically induced - wakefield 
processes. Conclusions are given in Section 4. 

2. THE STOCHASTIC WAKEFIELD SIMULATION 

Here we describe the 3D runs from our first wakefield ex- 
periments, full details of this as well as other set ups, tests and 
more results can be found in lFrederiksenI (|2Q08|) . 




Fig. 1 . — Illustration of our computational setup; GRB photons (left) move 
through the ambient CBM plasma. As the intensity of the light pulse in- 
creases, the weights (size) of the photons vary, rather than the number. Color 
designates that photons carry different energies (through h>) as well. 



A synthetic thermal GRB is delivered to a quiescent CBM 
plasma by adding photons (particles) to the computational do- 
main on a volume boundary, as illustrated in figure Fig. [T] 
according to a prescribed light curve and spectrum. 

The lig ht curve is ins pired by B ATSE data fitted to a F RED 
function dKocevski et a l. 2003; Rvde & Svenssonll2002h with 
parameters chosen for BATSE trigger 3891, r = 1.26 and 
d = 2.61. The spectral window coincides with the BATSE 
window, G [33.6keV, 3.36MeV] - after correction for an 
estimated redshift of z ~ 1 .68. The spectrum is assumed to be 
Planckian: N^(u,t) = {huy^F{u,t) = {hvy^By{T{t),t), for the 
photon spectral radiance. We use the Stefan-Boltzmann law 
to model the temperature, T{t) oc FitY^^. Our model FRED 
has duration tgrb = 200, corresponding to about 1% of T50 for 
trigger 3891. 

We evolve our synthetic GRB by adding a constant number 
of photons per cell carrying continuous weights prescribed by 
the burst functions. Photon energies are sampled from MC 
integration of the (time-dependent) Planckian. The photon 
spectral radiance is then 

F(iy t) ^^^""'^^ ^^^""^ 
n^iy.^t) oc ^ — (X Wi(u^t) => n^iu^t) oc w/(0 , 

^ /=i /=i 

using our assumptions about the thermal GRB. The last sum, 
split into assigning time- variable weights and frequency vari- 
ability to individual particles, gives a high degree of flexibility 
compared to conventional PIC codes. 

The CBM plasma is assumed to be a uniform density 
medium, consisting of a hydrogenic, fully ionized, moder- 
ate temperature plasma with Vth^e ^ 0.1c. We have adopted 
conventional PIC scaling, setting all relevant natural con- 
stants equal to unity, except the ion-to-electron mass ra- 
tio, which is Mi/ Me = 256. The initial number of parti- 
cles is 40/cell/species, and the grid resolution 100x20x4000; 
this translates to {Lx,Ly,Lz} ^ { 12(5^x2^^x500(5^} for our 



choices of length units and plasma density. The burst duration 
in light travel length is Lgrb ^ 64^^. A relatively flat volume 
aspect ratio is a compromise to ensure skin-depth resolution at 
the lowest possible grid size to give any 3D structure a chance 
of growing marginally without spending excess computing 
time. Boundary conditions are for the lower (upper) bound- 
aries: plasma particles are specularly reflected (thermalized) 
and outgoing photons are removed. Fields are damped in a 
layer AL^ ~ 0. 1 on the boundaries to avoid spurious wave 
interference. 

Compton scattering is conducted through MC binary inter- 
action, using the full Klein-Nishina expression. Photons and 
electrons are selected pairwise at random for scattering ac- 
cording to their weights. To achieve detailed balance, scatter- 
ing is done by splitting the particles into a new part carr ying 
the C ompton scattered fraction of the old part (Haugb oelld 
l2005h . and an unscattered remaining part. We may write 

^(-f,e),scatt = ^{-f,e),old-^{-f,e),new 5 

where the weights of scattered electrons and photons are equal 
pairwise. The particle scatter fraction per time step then sat- 
isfies 

where o^(Pe, P7) denotes the microphysical Compton scatter- 
ing matrix for photons and electrons in a given cell. For a 
scattering fraction of 10~^, 10^"^ physical particles are scat- 
tered away from a 'mother' particle of size 10^^, for example; 
this illustrates the flexible and inherent continuous nature of 
the PhotonPlasma code. 

3. RESULTS AND DISCUSSION 
3.1. Wakefield- Plasma Gradient Effects and Transients 

As the pulse interacts with the quiescent plasma, electrons 
are strongly forward Compton scattered. The plasma reacts 
electrostatically to restore the displacement by pulling the 
bulk of the electrons backwards anti-parallel to the burst prop- 
agation direction, while this is turn back-reacts on the pho- 
tons. The effects can be seen as the 'dip' from z = 1500, and 
to the left, in figures [2l and [3l where scatter plots in momen- 
tum sub-space {z, |Pz| } are shown for electrons and photons, 
respectively. The scatter plots are from runs with twice the 
resolution and half the density of those for the main simula- 
tion (Section 2 and figlU). The well resolved and longer skin- 
depth reveals strong electrostatic beating in the plasma paral- 
lel momentum. These beats last of order rt ^ Se/c ^ I0~^s in 
the CBM rest frame, and particles caught in the wake reach 
Lorentz factors upward of 7^ ~ 30. 

Forwardly accelerated plasma - above the yellow line in 
Fig. [2]- relaxes to \pz,e\ ^ ^rUeC in a pace proportional to the 
burst duration, in order to counter-act the bulk backwards flow 
(the 'dip') initiated by the photons. The plasma undergoes vi- 
olent accelerations, to maintain charge-neutrality. Burst pho- 
tons keep the electron population 'inflated' as long as pho- 
ton free energy and anisotropy is available in the plasma; 
far downstream the 'inflated' non-thermal population lasts 
through the duration of our simulations. 

The violent electrostatic acceleration of plasma electrons, 
due to the sharp gradient in photon pressure, will produce 
bremsstrahlung and synchrotron photons. In subsequent up- 
scattering these photons will eventually also affect the spec- 
trum (Barbi ellini et al .1120 06). More importantly, millisecond 
variability could - in our assumed environment - arise di- 
rectly from the enhanced density variations, scattering and 
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Fig. 2. — Electron Scatter Plot: Color designates weight in log-scale - 
blue particles have weights w^iue ~ 10~^wo, where wq is the initial weight. 
The yellow line coincides with zero bulk z-momentum, p^^^ =0. For clarity 
photons (black) are offset vertically by pz^^ = +6. 
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Fig. 3 . — Photon Scatter Plot: same run as Fig. [21 but for photons (in color). 
The initial burst (moving right) is visible from r ~ 1500 to f ~ 1300. Pink 
arrows: local photon upscattering. Green lines: energy threshold for 7+7 ^ 
e'^+e~ for photons pairwise oppositely outside (away from pz = 0) these lines. 
Electrons (black) are offset by pz,e = +6. Here, p^,^ = [hu]z/meC^ . 

possibly pair production (leading to higher densities and more 
scattering). Such a variability, which is comparable to the 
plasma frequency in a medium with density no ~ 10~^cm 
~^ - lO^cm"^, could therefore be self-f ceding and grow in 
strength and duration. We may be observing a weak seed to 
such spike in the spectrum at high energy, marked by the two 
pink arrows in Fig. [3l Such photons spikes could facilitate 
pair production, and modify the GRB spectrum as the burst 
traverses the CBM. 

3.2. Magnetic Field Generation 

Most intriguingly, we observe growth of a strong, large- 
scale magnetic field in the downstream wakefield gradient re- 
gion. The scenario is captured in Fig. |4l A cross section of 
the 3D volume is shown for electrons, photons, ions and the 
transverse magnetic field, (top-to-bottom), for the 3D run 
described in Section 2. Barely visible in the electron popula- 
tion is the same electrostatic spiky structure that shows in Fig. 
[21 but somewhat weaker due comparatively higher density and 
lower resolution. 

From the ion density and electron density panels it is seen 



that both species undergo instability to form confluent fila- 
ments in the plasma (although the inertia difference plays a 
role). The filamentation build-up is initiated by fluctuating 
EM fields in C BM plasma ahead of the burst. As described 
in Section 13.11 the plasma forms a two-stream anisotropy in 
reaction to the burst forcing. The fastest growing instability in 
this case is the two- stream filamentation mode - tantamount 
to the Weibel instability dWeibell [T959i) - which produces 
quasi- static magnetic fields. The effect is self-sustaining as 
long as the photon momentum anisotropy is sufficiently high. 

Furthermore, ponderomotive particle forcing, (x 
-ql(mujy^\/E^, is likely to be active in the strongly oscil- 
latory wakefield plasma. This drift force is independent of the 
charge sign; both electrons and ions drift in the same man- 
ner. Such effects has been rep orted by other authors ( Hos hind 
l2008l:lLiang & Noguchil2QQ7h to being responsible for the ac- 
celeration of particles due to forcing by Poynting flux pulses. 
The coherent nature of the plasma forcing by these authors, 
contrasts our stochastic forcing. 

Determining how the magnetic field will decay is a chal- 
lenge yet to be met; modeling a larger plasma volume is 
necessary to determine the ultimate temporal and spatial de- 
velopment. Resorting to lower-dimensional models will not 
solve the problem; the instability is inherently 3D in na- 
ture and cannot develop realistically in lower dimensions (cf . 
iFrederiksen et al.|[2QQ4|) . We also observe in Fig. |4l that fil- 
aments, current structures and field scales are limited by the 
computational volume - i.e. they are 'boxed-in'. Nonetheless, 
they survive relatively far behind the final tail of the photon 
wakefield (the GRB). Both electrons and ions are responsi- 
ble for the density filaments, the current densities are borne 
almost entirely by the electrons - again - owing to the iner- 
tial difference. We may therefore expect that as long as the 
photon anisotropy is "large", the field structures will survive. 

We proceed to speculate that such magnetized wakefields 
volumes could provide 'magnetic walls' against which trail- 
ing ejecta could gradually produce a build- up of shocks 
and particle acceleration, as the ejecta traverses the pre- 
conditioned circumburst medium. 

4. CONCLUSIONS 

We emphasize the importance of bridging the gap be- 
tween the (statistical) sub-Debye domain and the macroscopic 
Vlasov-/super-Debye domain in astrophysical plasma model- 
ing, opening the possibility of a detailed kinetic photon de- 
scription, and allowing the modeling of processes such as 
neutron transport and decay, neutrino streaming and transport, 
pair production etc. Frequencies for wave-particle collisions 
and for particle-particle collisions (e.g. Compton) are compa- 
rable in magnitude for the physical conditions we have con- 
sidered in this Letter. In effect 

VcoW CX LOpe = C~^5e , l^hm OC ncaj ^ i^coll ~ ^bin , 

where n is the photon number density, ar the Thompson cross 
section, / the burst light duration. Subscripts "coll" and "bin" 
denote "collective" (wave-particle), and "binary" (particle- 
particle) - respectively. 

Our main findings from these first simplified GRB -CBM 
wakefield interaction simulations, which we have conducted 
with the new PhotonPlasma hybrid scheme described 
briefly in Section 2, are: 

1 . Strong photon density gradients plays an important role in 
the dynamics of the CBM while being traversed by a GRB. 
We therefore lend support from our studies to suggestions set 
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Fig. 4. — Current structures for electrons (panel 1), photon pulse (panel 2), ions (panel 3) and (transverse, out-of -plane) magnetic field structure (panel 4) in 

lO'^Ar. The GRB photon pulse has light-length Lc = c/tgrb = 200 ^ 646 e. The computational 



our 3D Wakefield simulation, at simulation time t = lOlOcUpl 
volume has dimensions {Lx. 



: {40, 8, 1600} or, in electron skin-depths, {Lx 



[l2Se,2Se,500Se 



forth by (iBarbiellini et aLll2QQ6h to study in more detail these 
stochastically induced wakefield processes. 

2. The GRB spectrum is locally affected by the same gra- 
dient density induced effects, and this leads to local high en- 
ergy fluctuations in the GRB spectrum and light curve on sub- 
millisecond-to-millisecond time scales. 

3. A transverse magnetic field, quasi-static in the CBM frame, 
is generated in the wakefield of a GRB traversing a quiescent 
plasma. The field scales are likely limited by our numerical 
capabilities to follow the growth. 

It has been argued th eoretically (iBeloborodovl 120021: 
[Thompso n & Madaull2Q0Qh that the CBM ahead of the trail- 
ing forward shock could be loaded with a relatively high den- 
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sity pair plasma component. From this work, and from ob- 
serving Fig. [3l we see that accounting for pair production 
(7 + 7 ^ ^"^ + is imperative. In future work we shall add 
such dynamics, to test in particular whether density fluctua- 
tions could grow and eventually lead to externally produced 
shocks as a consequence of 7-ray preconditioning of the cir- 
cumburst medium. 
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